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We study the impurity effect in the time reversal symmetry (T) breaking superconductor 
based on the Bogoliubov-de Gennes (BdG) equations. In T- violating superconductors, spon- 
taneous currents are induced around the impurity. The current patterns around the impurity 
reflect the structures of the Cooper pairs. We investigate impurity problem numerically for 
two kinds of T violating superconductors (p x + ip y and d + is) and investigate the currents 
around the impurity. We also study the effects of the magnetic impurity in p-wave (p x + ip y ) 
superconductor, especially in view of the zero-energy crossing of energy levels related to the 
phase transition of the ground state. 
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§1. Introduction 

One of the current important issues of the superconductivity is the occurrence of the time-reversal 
symmetry (T)-violating state© It has been investigated in connection with the local behavior of the 
order parameter at the interface of the Josephson junctions and the surface of the material. Another 
context of T violating state is the unconventional superconductor where time reversal symmetry 
is violated throughout the whole materials. Such examples are seen in the heavy fermion super- 
conductors UPt3, Ui-ajTh^Beia (0.017 < x < 0.045)B and the layered perovskite superconductor 
Sr2Ru04.& In general the violation of time-reversal invariance is a property of the magnetism 
and indeed, it is found that the T-violating superconductor shows various magnetic propertiesBi) 
In such systems, inhomogeneities where the Meissner screening effects are insufficient allow the 
appeared spontaneous supercurrent. This happens, for example, at the surface of samples and also 
on the domain wall between the degenerate states of the superconductor. The defects of the crystal 
lattice, in particular impurities, can also lead to the appearance of the magnetic properties of the 
T-violating state, and it is known that the spontaneous suppercurrents are induced in the vicinity 
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of an impurity in such a systems. Such an induced current around the impurities can generate 
the local magnetic field in the superconducting state. The local magnetic field can be detected 
by a high resolution probe like spin polarized muonJ'i and indeed has lead to the observation of 
the intrinsic magnetism in the superconducting state of such T-violating superconducting systems. 
There are some studies about the induced current around the impurity in T-violating supercon- 
ductor .Bi'l 5 © We have already investigated a microscopic mechanism of such an induced current 
around the non-magnetic impurity with the Cooper pair of the type p x ± ip y , which is a candidate 
of the Cooper pair symmetry of the S^RuOjlU and has angular momentum L z = ±1. We see that 
the circular current generated around the impurities reflects the angular momentum of the Cooper 
pair. But still there are only few studies on how the induced currents around the impurity depend 
on the type of the Cooper pairs. So it is interesting to study the influence of the the structures of 
the Cooper pairs to the patterns of the induced currents around the impurity. Here in this paper 
we investigate the induced currents around a non-magnetic impurity with different types of Cooper 
pair, especially we pay attention to the angular momentum of the Cooper pair. Basically we can 
make a distinction between two kinds of T-violating superconducting states whose Cooper pair 
have an intrinsic angular momentum like p x ± ip y , d zx ± \d zy (L z = ±1), d x 2_ y 2 ± id xy (L z = ±2) 
and those that have no angular momentum like d + is. Such a difference of the structure of the 
Cooper pairs will affect the pattern of the induced currents around the impurity. We also study 
the effects of the magnetic impurity in p x + ip y superconductor and investigate the differences from 
the non-magnetic impurity. Especially we pay attention to the bound state and the ground state 
transition as well as the induced current. 

The organization of this paper is as follow. In §2, we solve the Bogoliubov-de Gennes equation 
numerically and investigate the induced current around the impurity. In §3, we study the magnetic 
impurity effect on the p x + ip y superconductor, §4 is the conclusion and discussion. 

§2. Pattern of the induced current around a non-magnetic impurity 

In order to investigate the impurity effect on the superconductor, we solve the Bogoliubov de- 
Gennes equation numerically. The formalism is given below. We start with the B.C.S Hamil- 
tonian on the 2D square lattice 



where t and U is the hopping between the nearest-neighbor (n.n) sites and the impurity potential 
respectively and \i is the chemical potential of the system. We set the non-magnetic impurity at 
the site 0. The real space order parameter Ajj is defined as 
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where Vij is the attractive pairing interaction between the site i and j. The Hamiltonian (2.1) can 
be diagonalized by means of the unitary transformation U 
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where Nl is the number of the lattice sites and e m (> 0) is the excitation energy of the Bogoliubov 
quasiparticles. The fermion operators 7 are related to the original electron operator aj i(T by 
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Then, the self-consistent equation for the gap Ajj (2.2) is given as 

\,j = ~ Y V hj\M* " N L +j,N L +rrMi,N L +m] 



(2.5) 



m=l 



We solve the equation (2.1) - (2.5) iteratively until the self consistent conditions are satisfied. In 
the T-violating superconductor, the spatial variation of the gap function induces the supercurrent. 
The current in the a directions at the site i is expressed as follow, 



Jia — W ^ (Mi+e a +N,rrM* i+N,m ~ i+e a +N,m,Mi+N,m) 



m=l 
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(2.6) 



First, we solve the problem in the case of rj± = p x ± ip y , which has an angular momentum L z = ±1 
if the systems has cylindrical symmetry. We take the nearest neighbor interaction for Vij and shift 
the chemical potential from half filling in order to stabilize the p-wave pairing. We assume the 
state 77+ = p x + \p y is realized globally in the absence of impurity. In the calculation we take the 
lattice size 27x27 with periodic boundary condition and we take the value Vij = 3t, U = lOOt, the 
chemical potential [i = —0.7t. 
The numerical result satisfies the relation, 

A'i,i+e x = — Aj+ eXj j, 

where e x and e y are unit vectors in the x and y direction respectively. This result shows that the 
odd parity state is realized and there are no even parity components within the numerical accuracy. 
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Each component of the p x and p y of the Cooper pairs at the site i are defined as, 

Ap>') = ^(A v ;+ ex - &i,i-ex) (2-8) 

A Py (i) = - (A M+ey - A M _ ey ) (2.9) 

Due to the spatial variation of the order parameter, r/_ = p x — ip y component is induced near the 
impurity. We determine the weight of the r]± for each site as, 

We depict the spatial dependence of the order parameter in Fig.l(a-c). Fig.(l-a) shows the r] + 
component which is dominant in the bulk and are suppressed around the impurity. In Fig.lb-c, 
we can see the admixture component which vanishes in the bulk and is induced around the 
impurity. The admixture component r/_ has phase winding 2, like — e 2ie r]_(i), reflecting the 

difference of the angular momentum between dominant and admixture components of the Cooper 
pairs. The BdG equations have bound state which is localized around the impurity. The energy 
level of the bound state is not zero even in the unitary limit, if the system has not the electron- hole 
symmetry (see Appendix 1). 

In Fig. 2 we show the induced circular current. The spontaneous circular current are induced 
around the impurity and the current pattern are the reflection of the intrinsic angular momentum 
of the Cooper pairs. The intrinsic angular momentum turns into a circular current around an 
impurity due to scattering. Within this numerical calculation on the lattice systems, it is difficult 
to observe the 'counter current' which we discussed in RefS) because it is small and extends beyond 
the lattice treated model. 

We also show another type of the T-violating superconductor. We take d xy + is, which is easily 
realized by the next-nearest neighbor interaction. In this case, the Cooper pair has no angular 
momentum (equal admixture of the orbital momentum L z = ±2) and is not a chiral state. We 
define the d-wave and (extended) s-wave components of the gap function as below, 

Arf(i) = — (Aj^-fex+ey + Aj^—ex— ey — Aj^—ex+ey — Aj^+ex— ey) 

A s (i) = — (Aj^+ex+ey + Aj i i_ ex _ e y + ^i : i~e-x.+ey + ^i,i+e-x—ey ) (2T1) 

In the calculation we take the next-nearest neighbor interaction Vi = V = 2t and the chemical 
potential jj, = —0.7t as in the p-wave case. In Fig.3(a-c) we show the self-consistently determined 
order parameters. The magnitude of the bulk d xy and s-wave component is estimated as A^ = 0.6t 
and A s = 0.22t within these parameters. Around the impurity the real parts of the s-wave order is 
generated. Not that also the s-wave component is influenced by the non-magnetic impurity, since 
it has a non-trivial phase dependence in the Brillouin zone. The real s-wave component is induced 
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Fig. 1. The self consistent gap function for the dominant r/+(i) component and the admixture r]-(i) components. 
Fig.(l-a) is the spatial variation of the »7+(i) component. Fig.(l-b) and (1-c) are the real and imaginary component 
of 7/-(i) respectively. 

along the diagonal direction and vanishes in the node directions of the d-wave gap, and change sign 
under 90-degree rotation as in previous studies, 0' which are reflections of the sign change of the 
d-wave component. The induced current around the impurity is shown in Fig. 4. In this case, the 
induced current is not circular but flows along the diagonal direction. The pattern of the current 
is understood from the Ginzburg-Landau theory. 
The Ginzburg-Landau free energy is given as, 

/ = / ^ME KIA.f + /3,|A J f + ^|nA J f) + 7l |A d | 2 |A s | 2 + i 72 (A*W + A (i 2 A* 2 s ) 

j=d,s 

+ K(U* x A s U x A* d - U* Y A s U Y A* d + c.c) + g s 5(r)\A s \ 2 + g d S(r)\ A d \% (2.12) 
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Fig. 2. Induced current around the impurity in p x + ip y superconductor. 

where g s and g d are the impurity potential for s and d-wave respectively and X and Y are directed 
to the diagonal (1,1) and (-1,1) respectively and II = — i(Vx, Vy). The current are expressed as, 

j(r) = (-2e)[(^ s A%n*A% + K d A* d U*A* d ) 
+ K{(A\W x A d + A* d U* x A s )e x + c.c} 

- K{{A* s H* Y A d + A* d nVA s )e y + c.c}] (2.13) 

where ex and ey are the unit vectors in the X— and Y— direction respectively. The parts which 
come from the mixing gradient term {K term ) is important for the induced current. The d-wave 
gap are most affected by the non-magnetic impurity, and the suppression of the gap occurs in the 
diagonal direction. So the current in the diagonal direction is most remarkable. Due to the sign 
change of the d-wave gap, the direction of the current is reversed between X— and Y— direction. 
If we consider the conservation of the current, the overall feature of the induced current should be 
given like Fig. 5, which we cannot see unfortunately them easily in the numerical data due to the 
their small value. But it is a natural feature because the superconducting state is admixture of the 
orbital angular momentum of Cooper pair, L z = 2 and L z = —2 ( and s-wave component L z = 0). 
The direction of the current circle is the appearance of the mixture of L z = 2 and L 2 = — 2 of the 
Cooper pair. 

The induced current is clearly distinguished from that of the chiral p x + \p y one. In chiral p x + ip y 
case, the induced current forms the magnetic moment around the impurity, but in the d xy + is 
case, the induced current is arranged like a quadrapole, which is just the reflection of the angular 
momentum of d xy component. We discuss the magnetic field generated by this type of Cooper pair 
(d + is) in Appendix 2. The d x i_ y 2 + is, state which are often discussed in the context of the high 
temperature superconductor, will give the same type of the induced current with 90 degree rotation 
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Fig. 3. The self-consistent gap function for the d xy + is superconductor. Fig.(3-a) and (3-b) shows the spatial 
variation of the d xy and s-wave components, respectively. Fig.(3-c) is the induced real component of the s-wave 
gap which is absence in the bulk. 

about the z axis. 

In this section we consider the superconductor which violates the T symmetry in the bulk. It is 
possible that the T-violating superconducting state in realized locally around the impurity© The 
induced current pattern will be the same as that of the corresponding T-violating bulk one. So if 
we can detect the magnetic field pattern around the impurity, we can distinguish the symmetry 
of the Cooper pair. But it is difficult to produce the locally T-violating superconducting state 
only by the non-magnetic or magnetic impurity without spin-orbit interaction. We can understand 
the reason in the GL treatment. Based on GL free energy (2.12) and we assume the bulk state is 
the real d-wave components only. The spatial variation of the order parameter induces the s-wave 
component near the impurity. If the impurity effect is not so strong, we expand the GL equation 
by the impurity potential. We can easily see the first order of the impurity potential, only the real 
part of the s-wave and d-wave component coupled to the impurity potential. So the phase of the 
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Fig. 4. Induced current around the impurity in d X y + is superconductor. 



Fig. 5. Schematic figure of the induced current around the impurity for d xy + is superconductor. 



induced s-wave is real with in the first order of the impurity potential and does not violate the T 
symmetry. 

§3. Magnetic impurity in p x +iPj/ superconductor 

In this section we discuss the magnetic impurity effects in chiral state, p^+ip,, superconductor. 
Here we limit our considerations to a classical spin, S » 1. In the singlet case we can take the 
direction of the impurity spin parallel to the z-axis by suitable unitary transformation. In the case 
of the triplet superconductor, we cannot perform the same operation due to the lack of the spin 
rotation symmetry. But we take the direction of the magnetic impurity in the z direction due to 
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the physical condition as we explains below. The spontaneous circular current will flow around the 
magnetic impurity and the magnetic field generated by the induced current will fix the direction of 
impurity spin to the z-axis. So we take the magnetic impurity as a classical Ising type and neglect 
the transversal components of the impurity spin. 

The Bogoliubov-de Gennes equation with the simple Ising-like magnetic impurity is given like 
below, 

h (r)ii i(T (r) + <rJS z 5(r)ii i(T (r) + J2a> I dr h a ^(r, r')u i _ (T /(r') = E i Ui <T (r) 
-ho(i>i_ CT (r) + o-JS z <5(r)u i _ (J (r) - / dr A* aa >(r,r )u ia '(r ) = E^^r), 

where J is the exchange energy and S z is the impurity spin. The BdG equation is decoupled by 
(«j(r), v j(r)) and (iij.(r), v^(r)) due to the absence of the transversal component of the imputiry 
spin. By the solution of (uj(r), v j_ (r)) with the eigenvalue Ei, we can construct the solution of the 
pair (it^(r), v-f(r)) with the eigenvalue — E; like (u*j p (r),it*j(r)) and this is an one to one correspon- 
dence. First, we neglect the spatial dependence of the gap function, and fix A(k) = Ao(k x + ik y ). 
The BdG equation (3.1) has two types of the solutions, one is the continuum state and the other 
is the bound state at the impurity. The bound state is degenerate and is explicitly given by, 

»B t (r>\ / JgfeW*-) W-^ iA (^)e«\ 

where No is the density of state at the Fermi level and c = 7rNoJS z and Io 2 = ~^( 1 ^ c a )^, and the 
functions fi(kpr) and f2(kF r ), where kp is the Fermi wave number, are approximately given for 
the distance r » 1/kp, as 



7T 7T V A 2 - g2 r 

/l( ^ r) "2\/^ COs( ^ r -4 )e 

(3.3) 

7T / 2 37T V A 2 -* 2 



e_ is an eigenvalue of the bound state and is given as 

A 



(3.4) 



We can see the mixing of the angular momentum of the wave functions in ( |3.2[ ). Here the above 
eigenvalue and the bound states are derived by the assumption of the electron-hole symmetry and 
the infinite band width of the conduction electron. We can easily see that for the electron-like state 
the particle wave function couples to the hole wave function of an angular momentum reduced by 
1. For the hole-like states it is just the opposite way around. This property is responsible for the 
fact that these states carry a circular current. The solution of the bound state with an eigenvalue 
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e + = is given as, («BJ.( r )) u Bt( r )) = i v *Bi( r ), ^*Bt( r )); as we said. The spectral density, 

A a (r = 0,uj) = a (0)5(uj — e) + Z^~> a (0)5(uj + e), at the impurity site due to the bound state is 
given like, Z^\ = Z^\ = 27rN " A f , and Z^\ = Z<+) t = 0. 
The current composed by the bound state is given below, 

■ b 4eA N |c| f 2 {k F rf 

36 ~ 77~1^ r ( ' 

7rm(l + c z J 2 ' 

In contrast to the non-magnetic impurity ,0) there is finite electron-like contribution of the current 
at the absolute zero temperature. This is due to the fact that the attractive nature of the exchange 
impurity potential (J > 0) for the down spin electrons. 

We can also calculate the contribution of the continuum state to the induced current as in the 
non-magnetic impurity case© In the unitary limit (c — ► oo), the current from the continuum is 
given by, 

eA , ,cos 2 (k F r - §) , - 2A o r , . , 

30 ~ -^h^ - 1 2 1 + e VF > ( 3 - 6 ) 

where vf is the Fermi velocity. The power-law decay 1/r 2 is because of the fact that the current is 
carried by the quasiparticle belonging to the continuum spectrum. 

Basically, the feature and the mechanism of the induced current is not different from the non- 
magnetic case and the contribution of the continuum spectrum is important. The direction of the 
current composed by the bound state is opposite to that of the continuum one. 
We can see the induced circular current also flows around the magnetic impurity and it does not 
contradict the assumption we first set up. 

Due to the spin degree of freedoms magnetic impurities cause another effect to the superconduct- 
ing state. In the s-wave superconductor, there is a zero temperature first-order phase transition 
depending on the magnitude of the impurity exchange potential J. The total spin of the electron 
change from spin unpolarized state, < s z >= 0, to the polarized one < s z >= — i, which was 
noticed by Sakurai.0) The physical picture of the phase transition is that the electron is captured 
by the magnetic imputiry due to the large enough exchange energy, J. So there is a competition 
between the pairing-condensation energy and the magnetic impurity interaction. The feature of 
the transition can be understood from the level crossing of the bound state. EE© We express the 
solution of the BdG equation of the pair (« n +(r), v n \ (r)) as 7 n and as j' m for the solution of the 
pair (w m j(r), u m f(r)). We set the n > for the eigenvalue E n > and n < for E n < 0. The 
difference of the number of up spin electrons and that of down spin electrons, M 2 is expressed as, 

M - = E 1 " E 1 + E 7 f n7n " E VWm (3-7) 

n<0 k n>0 m>0 

where the summation J2k runs over the wave number space. The spin of the ground state is given by, 
J2n<o 1 — J2k 1; an d changes its value when the number of the negative energy states are decreased 



10 



or increased, namely it is the level crossing across the Fermi level. We can see the energy level of 



the bound state for the p x + ip y superconductor is €b = — Ao/Vl + c 2 . It does not cross the Fermi 
level however large the exchange energy of the impurity is. So it can be concluded that a phase 
transition, as found in an s-wave superconductor, is absent in this case. This conclusion is also 
applicable for the d-wave superconductor and was already pointed out by Salkola© However, we 
used various assumptions like electron-hole symmetry and impurity s-wave scattering. These are 
rather ideal conditions for the actual magnetic impurity systems. If we give up these restrictions, 
the situation is changed. For example, if we allow the electron-hole asymmetry like N(£) = No+N'£, 
where N(£) is the density of state and N' is the introduced phenomenological parameter for the 
electron-hole asymmetry. Then the level of the bound state is given as, 

(1 - N'DJS Z ) a . , 

e± = ± . V z> =A (3.8) 



(ttN JS z ) 2 + (1 - 2N'DJS Z ) 2 

where D is the half of the band width. We can see that the level crossing can occur at the critical 
strength, J C S Z = 

We can also see that p-wave scattering at the impurity can also cause the phase transition to the 
spin polarized state. For example, we take the magnetic impurity potential as, 

J(k, k') = J s a z + J p (e^-^ + e-'^^Vz (3-9) 

where J s and J p is the s-wave and p-wave scattering term, respectively, and (f> and 4> is the angle 
of k and k' respectively. In this case the energy levels of the bound states solution for the pair 
(u n i(r),v n i(r)) are given like below (see Appendix3), 

(1 - cd)A 

£1 = ^(c + d^ + d-cd)* (3 ' 10a) 

• = "TTTP (3 ' 10b) 

where c = 7rNoJ s S and d = 7rNoJ p S. We can see from (3.10a) that the impurity spin can be 
compensated when c > 1/d. Another bound state is also accompanied by introducing the p-wave 
scattering (3.10b). We take another impurity potential like below, 

J(k, k') = U s + J p (eW-*'> + e-'^^Vz (3.11) 

where U is the s-wave scalar potential. In this case, the levels of the bound state is given like, 

1 — cd 

ei = - A (3.12a) 

y/{c + d) 2 + (1 - cd)^ 1 ; 

£a = 'A (3 - 12b) 

63 = ± v^WWTW Ao (3 - 12c) 
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where c = 7rNoU s and in (3.12c) we take plus sign for (c > d) and minus sign for (c < d). The form 
of (3.12a) can also be possible for the level crossing. In the unitary limit for the scalar potential 
c — > oo, the level crossing occurs at infinitesimal p-wave magnetic potential. The bound state energy 
level with the unitary limit in s-wave scattering, c — > oo, is idAo/Vl + d 2 and does not coincide 
with the zero energy level. The above scenario is applicable to the d-wave superconductor with 
d-wave scattering from the impurity. The p-wave and d-wave scattering is given naturally from the 
nearest neighbor exchange interaction of the impurity like Jo J2t So a o+r, EE) where summation 
r is taken at the nearest neighbors of the impurity site and Sq is the impurity spin and a the 
electron spin. Basically, the impurity resonance peak at the zero energy level in unitary limit splits 
by many factors like electron-hole asymmetry and p-wave scattering. The resonance peak splitting 
of the DOS in Tsuchiura© may be related to the anisotropic scattering from the impurity. The 
anisotropic scattering from the impurity is not negligible in actual transition metal systems like 
the cuprate superconductors. So actually there is also a first order phase transition of the ground 
state in this p-wave case as in the conventional s-wave superconductor. But the critical strength of 
the impurity exchange potential J c is not determined only by the competition of the condensation 
energy and the exchange energy as in the s-wave case but various factors should be considered. 
We carry out the numerical calculation of the BdG equation as in the non-magnetic impurity. In 
Fig. 6 we show the spatial dependence of the gap function. Basic feature is not different from that of 
the non-magnetic impurity case. The magnetic impurity induces the spatial dependence of the spin 
density. In Fig. 7 we show the spin density around the magnetic impurity. The peak < s(r = 0) > 
is cut off in order to illustrate the fine details. In Fig.(7-a) the impurity potential JS = 5t and is 
not the electron spin-unpolarized state (< s z >= 0). On the other hand, (7-b) is the spin-polarized 
state (< s z >= — |) with JS = lOi. There is a compensating spin-density cloud of opposite sign in 
the neighborhood of the magnetic moment in order to screen the spin polarization as in the d-wave 
case.0) On the other hand, for the spin-polarized state, the compensating spin-density cloud is 
much smaller than that of unpolarized one due to the imbalance of the spin up and down electron. 
Just at the impurity site, r = 0, the bound state contribution to the spin density < s(r = 0) > 
changes at the transition point from the up spin to the down spin. So there exists the discontinuity 
of < s(r = 0) > at the transition. But the main contribution of the < s(r = 0) > comes from the 
continuum part of the spectrum. 

§4. Conclusion and Discussion 

In this paper we have investigated the non-magnetic impurity effect by solving the BdG equation 
numerically for chiral p x + \p y and the non-chiral but T- violating d xy + is superconductor in §2. 
We see induced currents of the p x + ip y and d xy + is superconductors. Both of the T- violating 
superconductor can generate the induced spontaneous currents around the impurities. The induced 
current can generate localized magnetic field in the superconductor. But the patterns of the induced 
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Fig. 6. The self-consistent gap function for the p x + ip y superconductor with the magnetic impurity. The calculation 
is done with the parameter /i = — 1.5t and JS = 5t on 27x27 lattice. Fig.(6-a) is dominant p x +ip y component and 
(6-b), (6-c) is the admixture real and imaginary p x — ip y component, respectively. 

currents are quite different reflecting the property of the Cooper pairs. In the chiral superconductor 
Px + ipy, circular currents around the impurity are induced spontaneously. On the other hand, in 
the d xy + is superconductor, the current pattern has quadrapolar shape reflecting the mixture of 
the angular momentum L z = 2 and L z = —2. 

In §3 we investigate the magnetic impurity in p x + ip y superconductor. In this case, the circular 
current is induced as in the non-magnetic impurity case. But the bound states also carry the 
current in contrast to the non-magnetic case. The level of the bound state is not different as 
in the non-magnetic case and the phase transition of the ground state is absent when the 'ideal' 
conditions, such as electron-hole symmetry and s-wave scattering from the impurity, are satisfied 
in the p-wave case. But if some conditions are violated, then there is a transition to the electron 
spin-polarized state. Especially, we pay attention to the anisotropic scattering from the impurity 
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Fig. 7. Spin density around the magnetic impurity. The strength of the magnetic impurity JS is 5t and 8t for (a), 
and (b) respectively. Total spin moment is (a) < s z >= and (b) < s z >= — \. 

which cause the phase transition to the spin polarized state. 

In this paper, we take a classical impurity spin and neglect the transversal components of the 
magnetic impurity spin. So we do not consider the Kondo effect is the usual sense. It is expected 
that in the T violating superconductor the induced current pins the impurity spin to the parallel 
direction of the generated magnetic field. This effect would weaken the Kondo effects in addition 
to the presence of the superconducting gap. But if the Zeeman energy is small, then the Kondo 
effect should be taken into account. It is interesting to investigate the competition of the Kondo 
effect and this Zeeman energy in T-violating superconductor. Other than this effect, in the triplet 
superconductor, there are many interesting effects in connection to the Kondo effect. Especially 
due to the lack of the spin symmetry of the superconducting state and the mixture of the finite 
angular momentum channels other than s-state (angular momentum m = channel), the Kondo 
effect in the triplet superconductor may have different features from the s-wave case. This problem 
is left for future studies. Also the spin-orbit interaction shows interesting effects. The spin-orbit 
interaction directly couples to the orbital part of the Cooper pair and give different effects for the 
p x + ipy and p x — ipy state. So the induced components, p x — ip y , around the impurity in p x + \p y 
superconductor will be affected by the spin-orbit interaction. Also this problem will be discussed 
elsewhere. 
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Appendix 1 

Electron-hole asymmetry shifts the the energy levels of the bound states. Here we concentrate 
the p x + ip y type of Cooper pairs. As for impurity potential, we take the contact type s-wave 
scattering like U5(r). The Fourier transformed BdG equation is written as, 

ikUka + AkVk-a +Tf^2 U ka = eU ka (4.1) 
k 

-CkVk-a + A* k Uka ~ T- V k-a = £V k -a (4.2) 

k 

where u k and v k are electron and hole part of the wave function, respectively, is quasi-particle 
dispersion and is gap function. We introduce the variable Io = (U/V") J2k u ka and Io' = 
(U/V) J2k v k-a, and the self-consistent equation for Io and Io' are given as 

r _U r (e + &)I - Afclp' _ T < ■ j,, 



V k a ~ V k e 2 "6 2 -|A fc | 2 " i0 V^e*-tf-\A k \i 



(4.3) 



where for final form we used the fact that the angular integral over the gap function A k vanishes. 
The £fc term in the numerator does not vanish due to the electron-hole asymmetry. The solutions of 
( [Q| ) are divided into two types (Io / and 1$ = ) and (I = and I ' ^ ). Here we introduce 
the electron-hole asymmetry phenomenologically like N(£) = No + N'£, where N(£) is the density 
of state. Then with neglecting the fc-dependence of the |A^| 2 in denominator of (|4.3| ), which are 
approximately valid in the p x + ip y case, in the unitary limit U — > oo, the energy level of the bound 
state ( e < |Afc| ) is given like, 

e ~ ± ' 1 (4.4) 

^/(Novr) 2 + (2N'D) 2 

where D is the half of the band-width. The energy level does not coincide to zero energy level but 
splits due to the electron-hole asymmetry. The above argument is also applicable to the d-wave 
superconductor. 
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Appendix 2 

In this Appendix we discuss the magnetic field generated around the impurity in the d + is type 
Cooper pair. Our discussion is based on the GL theory. In the case of the coexistence of the d and 
s wave with T-violating superconductor, the GL free energy in the weak coupling approximation 
is written as 

/ = -2a s |A s | 2 - a.lA,! 2 + /?(|A S | 4 + jj|A d | 4 + 2|A s | 2 |A rf | 2 + ^[A s * 2 A d 2 + A S 2 A/ 2 )] 
+K[2|nA/| + |IIA/| + (n/A s n x A/ - U y *A s U y A d * + c.c)] + 5 (|A S | 2 + \A d \ 2 )5(r) 

(4.5) 

where a s = ln^p- and a d = ln-^ and we set the same strength of the impurity potential for d and 

s wave order parameter. 

The derived GL equations are, 

-2aS + ^(\S\ 2 + \D\ 2 )S + ^D 2 S* + 2U* 2 S + (IT^ - U 2 y )D + gS5(r) = (4.6) 

. 

3M 1 1 1 1 ' 3 

j = [2S*US + D*UD + (SU X D* + DU x S*)x - (SU y D* + DU y S*)y + c.c] (4.8) 



—D + -( 7 |£>| 2 + 2\S\ 2 )D + -S 2 D* + U* 2 D + (II 2 s - U 2 y )S + gD5(r) = (4.7) 



here a = ^ and II = — iV. D and S is the order parameters for d- ans s-wave, respectively, which 



are normalized by Ao = y and the length scale is also normalized by y ^ . For the realization 
of the bulk T-violating state we impose the condition 1 > a > |. We expand S and D in first 
order of impurity potential, g, like S = i(So + Si) and D = (Dq + Di), where So and Dq are the 
values of the homogeneous case and are given as So 2 = f (3a — 2) and Dq 2 = 3(1 — a). The Fourier 
transformed form of Si and Di are rather complex and are given like 

, , 2 V-6 + 9 a c/ (16 a 2 + 3 q 2 (-4 + q 2 ) - 2 a (8 - 9 q 2 ) - q 2 (4 - 2 a + q 2 ) cos(4 0)) 
o (r/. o) = 



S (q,cf>) = 
D'(q,4>) = 



C{qA) 

4 V3 - 3 a g cos(2 <ft) (16 + 48 a 2 - 16 q 2 + 3 q 4 + 8 a (-7 + 3 q 2 ) - q 4 cos(4 <p)) 

Cfe0) 

8V3-3ac/ (-16a + 24a 2 -10q 2 + 21ag 2 + 3q 4 -g 2 (2 - a + q 2 ) cos(4 0)) 

C(^ 



4V-6 + 90C/ cos(2 0) f-2 (8 + 16 a 2 - 6 q 2 + q 4 + 6 a (-4 + q 2 )) + g 4 cos(2 0) 2 ) 

D <«■ *> 

where we set (q x ,Qy) = (<?cos(</>), qsm(4>)) and decompose like Si = S' + iS 1 " and Di = D' + \D ". 
The denominator factor, C(q,<f>), is given as, 

C(<7, <p) = 512a - 1280a 2 + 768a 3 + 320<? 2 - 896a<? 2 + 480aV - 48<? 4 - 48aq 4 - 19<? 6 

+ Aq 2 (l6 + 24a 2 + 4g 2 + 3g 4 + 4a (-8 + q 2 )) cos(4<£) - q 6 cos(80) (4.10) 
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The current in the first order of g is given as 

j(r) = (-4S d x S"(r) - 2D d x D" '(r) + 2S d x D' (r) - 2AAS'(r))x (4.11) 

+ (-4So^5"(r) - 2 J D 5 2/ £»"(r) + 2S <9 J/ £>>) - 2D 8 y S' (r))y (4.12) 

The magnetic field can be calculated by the Biot-Savart formula like, 

q x j(q) 



B(g) 



9 2 



q x q y (S D\q)-D S'(q)) 



q 2 



(4.13) 



We can easily see the sign of the magnetic field changes by 90 degree rotation in (q x , q y ) space from 
(4.9) and (4.13), so the feature of the generated magnetic field is the quadrapole like. The feature 
is consistent to our numerical calculations in §2. 

Appendix 3 

The energy levels of the bound states with p-wave scattering from the impurity are determined 
as below. We take the same form of the magnetic impurity potential as (3.9). Then as in the 
Appendix 1 from the Fourier transformed form of the BdG equation we get 



, M (ejv + e^ao + Iie^ + Iie-^ + A^Jo + Jie^ + Jie-^; 
u tN (k) = £2n _ ^ _ Aq2 (4.14a) 



m A* fc (I + Iieg + 1'ie-^) + (eg - &)(J + Jie^ + JV 1 *) 
viMk) = £2n _ ^ _ Aq2 , (4.14b) 

here we define as I = J S S J2k' u ],N(k'), h = J P S J2k' ' u T,A r (^') e ~ 1 * > = J pS J2k' ^LaK^'K^ > 
Jo = JsS^fe' v l,N(k'), Ji = JpS^fc' v l,N(k')e~ lc ^ and j[ = J P S vi,n (k')e 1 ^ , where ^' is the 
angle of k' . Assuming |ejv| < Ao, we get the bound energies from the linear equations of Io,Ii,l'i,Jo, 
Ji and j'i. The energy levels of the bound states with the impurity potential (3.11) are given by 
replacing Jo to — Jo in the above expression (4.14). 
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